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ABSTRACT 

Cosmic rays that escape their acceleration site interact with the ambient medium and 
produce gamma rays as the result of inelastic proton-proton collisions. The detection of 
such diffuse emission may reveal the presence of an accelerator of cosmic rays, and also 
constrain the cosmic ray diffusion coefficient in its vicinity. Preliminary results in this 
direction have been obtained in the last years from studies of the gamma-ray emission 
from molecular clouds located in the vicinity of supernova remnants, which are the 
prime candidate for cosmic ray production. Hints have been found for a significant 
suppression of the diffusion coefficient with respect to the average one in the Galaxy. 
However, most of these studies rely on the assumption of isotropic diffusion, which 
may not be very well justified. Here, we extend this study to the case in which cosmic 
rays that escape an accelerator diffuse preferentially along the magnetic field lines. As 
a first approximation, we further assume that particles are strongly magnetized and 
that their transport perpendicular to the magnetic field is mainly due to the wandering 
of the field lines. The resulting spatial distribution of runaway cosmic rays around the 
accelerator is, in this case, strongly anisotropic. An application of the model to the 
case of the supernova remnant W28 demonstrates how the estimates of the diffusion 
coefficient from gamma-ray observations strongly depend on the assumptions made 
on the isotropy (or anisotropy) of diffusion. For higher levels of anisotropy of the 
diffusion, larger values of the diffusion coefficient are found to provide a good fit to 
data. Thus, detailed models for the propagation of cosmic rays are needed in order to 
interpret in a correct way the gamma-ray observations. 

Key words: cosmic rays - gamma rays - ISM: supernova remnants. 



1 INTRODUCTION 

Galactic Cosmic Rays (CRs) are mainly constituted by rela- 
tivistic protons and are believed to be accelerated at Super- 
Nova Remnant (S NR) shocks via first order Fermi mecha- 
nism |Hillasil2005h . Though very popular, this scenario still 
needs to be conclusively proven by observations. 

If CRs are indeed accelerated at SNRs, these objects 
must be gamma-ray sources. This is because the CRs ac- 
celerated at the shock undergo inelastic proton-proton in- 
teractions with the ambient medium and pro duce neutral 
pions which in turn decay in to gamma rays l|Drurv et al.l 
ll994l ; lNaito fc Takaharalll994f ). Several SNRs have been de- 
tecte d in gamma ray s at both TeV ( e.g. iHinton fc Hofmannl 
120091 ) and GeV (e.g. iGiordandlioill ) energies, in agreement 
with such expectations. However, it is often difficult to de- 
termine whether the origin of the gamma-ray emission is 
hadronic, and thus related to the acceleration of CRs, or 
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due to leptonic mechanisms such as inverse Compton scat- 
tering. For this reason, multi-wavelength studies of SNRs 
have been extensively carried out in an attempt to solve 
this degeneracy. Though for some individual SNRs it has 
been possible to ascribe the gamma-ray emission exclusively 



and quite confidently to hadronic (e.g. lAcciari et al 
Morlino fc Capriolill2012h or leptonic fe.g. lAbdo et a. 



2011; 



2011 



Ellison et alH2010l ) processes, in other cases this ambiguity 



remains a problem. 

An alternative way to reveal the presence of a CR 
source is by searching for the r adiation produced by CRs 
that escape the acceleration site jXh aronian & Atoyan ji996l ; 
Gabici fc Aharonianl 120071 ; iRodriguez Marrero et al.l 120081 ; 



Gabici et al.ll2009f) . At some stage of the dynamical evolu- 



tion of the SNR, CRs are expected to leave the shock region 
and escape into the interstellar medium. The details of the 
esca pe mechanism are still not very well understood (see 
e.g. iGabicil l201ll . and references therein) , but it is gener- 
ally believed that the ability of a SNR in confining parti- 
cles decreases gradually with the shock speed, with higher 
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energy particles leaving the shock earlier than low energy 
ones. Once escaped, CRs diffuse away from the SNR and 
produce gamma-rays in interactions with the ambient gas. 
To date, some possible evidence for particle escape from 
SNRs has been pushed forward by the observation of dif- 
fuse gamma-ra y emission from the vicinity of the shell of 
the SNRs W28 dAharonian et alJliooj; iGiuliani et al.ll2ofol ; 



lAbdo i et al]|2010h and W44 jUchivama et alj|2012j) . In both 
cases, the emission is clearly located outside of the shell and 
it is spatially coincident with the location of massive Molec- 
ular Clouds (MCs). This would favor a scenario in which 
the MCs are illuminated by the runaway CRs and, being 
very massive, become p rominent gamma-r ay sources (for 
the case of W28 see e.g. iGabici et af]|201Cl and discussion 
in Section 

Besides revealing the presence of a CR source, the 
gamma-ray emission from runaway CRs can also be used 
to constrain the particles' diffusion coefficient in the re- 
gion surrounding the accelerator. This is very important 
for several reasons: first of all, a theoretical determina- 
tion of the diffusion coeffic i ent is a very complex task (see 
e.g. lYan fc Lazarianl I20041 . 120081 ') and observational con- 
straints are needed in order to guide and constrain mod- 
els. In addition to that, the diffusion of CRs is believed 
to be a non-linear process in which CRs themselves gen- 
erate vi a streaming instability t he turbulence they scatter 
off (e.g. iKulsrud fc PearceHl969T l. This is particularly rele- 
vant close to CR sources, where the CR density is expected 
to be very high, and possibly sufficient to suppress signifi- 
cantly the diffusion coefficient through st reaming instability 
|Ptuskin et al.ll2008l ; iMalkov et al.ll2012h . Thus, an empiri- 
cal determination of the diffusion coefficient can reveal pre- 
cious information on the ways in which particles and waves 
interact in astrophysical plasmas. 

Most of the studies aimed at the determination of 
the CR diffusion coefficient from gamma-ray observations 
rely on the assumption of isotropic diffusion ([Torres et al.l 
20081; iFuiita et al.ll2009l; iLi fc C hen 20101; IGabici et al.ll2010l ; 



Ohira et al.ll201ll ; lYan et al.ll2012h . The common rationale of 



these approaches can be summarized as follows: if SNRs are 
the sources of CRs, they have to convert a fraction n ~ 10% 
of their explosion energy Esn = 10 51 -Bsierg into accel- 
erated particles. If the diffusion of CRs proceeds isotrop- 
ically, after a time t from escape CRs of a given energy 
E are distributed roughly homogeneously within a distance 
Rd(E) « y/6D(E) x t from the SNR. Here, D(E) is the en- 
ergy dependent diffusion coefficient of CRs. Gamma-ray ob- 
servations of MCs located in the vicinity of W28 or W44 tells 
us which is the CR density ncn(E) needed to explain the ob- 
served emission. According to what said above, such density 
has to be of the order of ticr(E) « f sp (E) tjEsn / R 3 , where 
the factor f ap (E) contains the information on the shape of 
the spectral energy distribution of escaping CRs. For aged 
SNRs such as W28 and W44, the time t after the escape 
can be identified with the SNR age t aS e, and thus an ex- 
pression for the diffusion coefficient can be obtained, and 
reads: D « (f sp riEsN /ncR.) 2 ^ 3 tZge- Since the values of all 
the physical quantities present on the right side of the equa- 
tion can be inferred from observations, the expression pro- 
vides a direct estimate of the diffusion coefficient. 

As an example, w e summarize here the results obtained 
bv lGabici etafl (|2010h in interpreting the gamma-ray emis- 



sion observed from the MCs located close to the SNR W28. 
They obtained a good fit to the gamma-ray spetrum mea- 
sured by H.E.S.S. at photon energies > 300 GeV by assum- 
ing a diffusion coefficient for w 3 TeV CRs of the order of 
D(3 TeV) r; 5 x 10 27 (ry/0.1) 2/:i cm 2 /s. The corresponding 
diffusion length Rd of these particles is of the order of 100 
pc. Also the broad band gamma-ray spectrum from GeV to 
TeV energies can be fitted by adjusting the energy depen- 
dence of the diffusion coefficient and the distances between 
the SNR and the clouds. The important point here is the fact 
that the estimated diffusion coefficient is more than one or- 
der of magnitude smaller than the one normally adopted 
to describe the propagation of CRs of energy > TeV in 
the galactic disk, which is w 10 29 cm 2 /s (|Strong et al.ll2007l ; 
ICastellina fc Donate! I2012T ). These results are very similar 
to the one s obtained by other authors by means of similar 
modeling jFuiita et aill2009l ; ILi fc Chenl l2010l : lOhira et ail 
l201ll ; lYan et al.ll2012T l and seem to point toward a drop of 
the diffusion close to the SNR W28. 

However, the validity of the assumption of isotropic dif- 
fusion of CRs needs to be discussed. In fact, if the intensity 
of the turbulent field SB on scales resonant with the Lar- 
mor radius of particles is significantly smaller than the mean 
large scale field Bo (i.e. if SB /Bo <C 1), then CR diffusion 
becomes anisotropic, with particl es diffusing prefer entially 
along the magnetic field lines (e.g. lCasse et al.ll2002r i. In the 
limiting (but still reasonable) case in which the perpendicu- 
lar diffusion coefficient can be set equal to zero, the transport 
of CRs across the mean field is mainly due to t he wandering 
of magnetic field lines (|jokipii fc Parked 1 19691 ). This is the 
situation that we investigate in this paper. 

To give a qualitative idea of the role that anisotropic 
diffusion can play in these study, let us consider an ideal- 
ized case in which particles that escape a SNR diffuse along 
a magnetic flux tube characterized by a very long coherence 
length (i.e. the magnetic flux tube is preserved for a long 
distance). In this case, after a time t particle will diffuse up 
to a distance R d w <J2D\\ x t along the tube (here D\\ is 
the parallel diffusion coefficient), while their transverse dis- 
tribution will be equal to the radius of the SNR shock at 
the time of their escape, R s h, which is of the order of w 1— 
10 pc. Thus, the enhanced CR density in the flux tube will 
be proportional to ucr oc (RdR^)' 1 instead of oc R~/ 3 as 
in the isotropic case. It is easy to see that the estimates of 
the diffusion coefficient based on the two opposite assump- 
tions of isotropic and one-dimensional diffusion will differ by 
a factor of « (R d /R s h) 4/3 , which can be much larger than 
an order of magnitude! Thus, it is of paramount importance 
to investigate how the interpretation of gamma-ray obser- 
vations depends on the assumptions made concerning CR 
diffusion. 

In Section [2] we develop a model for CR propagation 
in which CRs are strongly magnetized and diffuse uniquely 
along the magnetic field lines. The wandering of the field 
lines is also taken into account, and a diffusion coefficient 
D m for the magnetic field lines that depends on the prop- 
erties of the turbulent field is defined (see e.g. iDuffv et al.l 
1995). In Section [3] the model is used to predict the spatial 
distribution of runaway CRs and their spectrum. Finally, we 
apply the model to the case of the SNR W28 in Section U 
A good fit to data is obtained, and the estimate of the par- 
allel diffusion coefficient is found to depend on the level of 



Anisotropic CR diffusion 3 



anisotropy of the diffusion. For higher level of anisotropy, i.e. 
smaller values of D m , larger values of the particles' diffusion 
coefficient are needed in order to fit data. A discussion of 
the results and of future perspectives in this line of research 
can be found in Section [5] 



2 COSMIC-RAY TRANSPORT IN THE 

PRESENCE OF MAGNETIC FIELD LINE 
WANDERING 

Consider a magnetic flux tube whose mean magnetic field 
Bo is assumed to lie along the z-axis, perpendicular to 
the (a;, y) plane. The wandering of magnetic field lines is 
due to long wavelength perturbations, i.e. perturbations on 
scales much larger than the particles' gyroradii, with root 
mean square amplitude SB. The condition for the valid- 
ity of quasi-linear theory is SB /Bo <C Lx/Ln, where L± 
and L|| are the field cohere nce lengths perpendicular an d 
parallel to Bo, respectively l|Kadomtsev fc Pogutsd [T979ft . 
According to quasi-linear theory, the field lines passing 
in the vicinity of (xo,yo) at z — are spread over a 
larger region as they reach a given z. The probability dis- 
tribution describing this spreading of field lines is gaus- 
sian and characterized by ((a; — xq) 2 ) = {(y — yoj 2 ) = 
2D m z, where brackets indicate an ense mble average and 
Dm, i s a diffusion coefficient for field lines (|Jokipii fc Par ken 
1969). For a broad band Fourier spectrum of the pertur- 
bation, the coherence lengths can be expressed as L±« = 
27r/Afcx,|| ~ 27r/fcx,||, with k± and fcii the characteristi c 
wave-vectors of the perturbation l|Achterberg fc Ballll 1994ft . 
Under these circ umstances, the diffusion coeffi ci ent is D m = 
(SB/ BpfLn/A (|Kadomtsev fc Pogutsd 1 19791 ; iDuffv et all 
1995). It is also possible to define the Lyapunov length 
\l = L\{5B I Bo)~ 2 1 Ln which describes the exponential 
separation of field lines w hose initial separation is smaller 
than L± (|lsichenkd ll99lT ). This can be interpreted as the 
length above which the flux tube is disrupted by field line 
divergence. For fiducial values of the parameters the length 
of t he flux tube is of the order of a few hundred parsecs (see 
e.g. iPtuskin et aIll200Sft . 

We are interested here in studying the propagation of 
CRs in the presence of magnetic field line wandering. Ener- 
getic particles diffuse along and across the magnetic field line 
as the result of resonant interactions with magnetic field per- 
turbations. Such perturbations are characterized by length 
scales of the order of the particles' Larmor radii. Such scales 
are much smaller than the ones responsible for field line 
wandering. According to quasi-linear theory, the ratio be- 
tween the parallel to perpendicular diffusion coefficient is 
Du/Dx = 1 + (A || ) 2 , where Ay is the particle's mean 
free path along the field line and r g is its gyroradius. In 
the interstellar medium i t is believed that An 3> r g which 
implies D± < D\\ (e.g. ICasse et all 120021 ). Thus, in the 
following we will neglect the diffusion of particles perpen- 
dicular to the field lines. In other words, a given particle 
remains attached to the same field line. Under this con- 
ditions, in a time interval At a particle diffuses along a 
given field line a distance ((Az) 2 ) — 2DuAt, but over such 
a distance Az along the z-axis the field line is displaced by 
an amount ((Aa;) 2 ) = 2D m Az . This leads to dGetmantsev 
1 19631 ; iRechester fc Rosenbluthl 1 19781 ; IChuvilgin fc Ptuskin 



Il993h : 



((Az) 2 ) oc D mS /D^Ki 



(1) 



which describes a sub-diffusive transport of particles per- 
pendicular to the mean magnetic field Bo- 

This behavior of energetic particles due to the combi- 
nation of the particle diffusion along the field lines and the 
random walk of the field lines themselves has been often 
referred to as compound diffusion, or anomalous diffusion. 
Models of compound diffusion have been developed and used 
in a great variety of context s, to study phenomena like the 
heat transport in Tokamak (|Rechester fc Rosenbluthlll978l ; 
Ilsichenkdll99lll. the propagation of e nergetic particle s in th e 
solar wind (|Jokipii fc Parked QUI; IZimbardo et~ai1 12006ft 
the co nfinement of CRs i n th e Galaxy jGetmantsev 
19631 ; iLingenfelter et all Il97ll ; IChuvilgin fc Ptuskin 



1993) and their acceler a tion at as t rophy si cal shocks 
jAchterberg fc Balll 1 19941 ; IDuffv et al l Il995l ; iKirk et all 
Il996ft . In this paper, we apply the formalism of compound 
diffusion to another context, which is the propagation 
of CRs in the vicinity of their sources, i.e. SNRs, after 
they escape the acceleration region. We will show that in 
this situation an accurate modeling is needed in order to 
interpret in a correct way the gamma-ray observations of 
MCs located in the vicinity of SNRs. 

In order to describe the compound diffusio n of CRs we 
adopt the mathematical formalism developed bv lWebb et al.l 
(2006) and we define Pyhw(x\z) as the probability to find a 
field line displaced by an amount Aa; after a step of length z 
along the direction of the umperturbed field Bo- From what 
said above, it follows that: 



r(Ax\z) = 



: exp 



(Aa 



4A, 



(2) 



which corresponds to a diffusive behavior of field lines (FLW 
stands for Field line Random Walk). A similar equation 
holds for the displacement Ay. This has to be combined 
with the probability P||(z|Ai) that a particle moves a dis- 
tance z along the field line in a time At — t—to- For diffusive 
transport of particles along the field we have: 



P\\{z\At) = 



y/4irD }l At 



exp 



4D„At 



(3) 



The probability for a particle to reach the position (x, y, z) at 
the time t, when its position at the time to was (a;o, yo, zo = 
0), is then the product of Pfrw with P«: 

P{Ax,Ay,z-At) = P||0|At) P FRW (Aa-|z) P FRW (Ay|z) 

(4) 

In order to m odel the escap e of C Rs from a SNRs we 
assume, following I Ptuskin et al.l (2008), that particles are 
injected in the flux tube in the xy-plane at z = 0, within 
a circular region whose radius is equal to the SNR shock 
radius R a h(E). Since CRs of different energy are expected 
to escape the remnant at different times, the radius of the 
injection reg i on is a n energy dependent quantity. Following 
iGabici et all (2009) we assume a power law scaling to con- 
nect the particle energy of the runaway CRs with the time 
after the supernova explosion: 



Em ax 



tSed 



(5) 
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where the implicit assumption has been made that the max- 
imum energy of CRs accelerated in a SNR Em ax is reached 
at the time tsed which marks the transition between the free 
expansion and the Sedov phases of the SNR evolution, and 
that CRs are gradually released in the interstellar medium 
from that time on. The Sedov phase is characterized by a 
scaling R s h oc t 2//5 which gives: 



Rsh{E e3C ) oc 



/ E eac \ 



5» 



(6) 



\Emax J 

which is what is assumed in the following. Other parameteri- 
zations of the escape time of CRs can be easily implemented. 

The spatial distribution of CRs can now be obtained by 
integrating the probability function given by Eq. [4] within 
the range Rq = yj (x 2 , + y%) < R 3h (t e3C (E)). To do so, it 
is convenient to adopt a cylindrical coordinate system and 
express the Field line Random Walk part of Eq.|4]as a func- 
tion of the quantities R — y/x 2 + y 2 , Ro = \J (x 2 + y 2 ) and 
cos(Aip) = (xxq + yyo)/(RRo) which leads to: 



f CR {R,z,t,E)=A 



E~ 



P«(z\t-t e 



dip 



d7? -Ro Pfrw{R, Ro, &<p\z) 



(7) 



where it has been assumed that the total spectrum of CRs 
released in the interstellar medium during the whole life of 
the SNR is a power law AE~ T with normalization: 



A = 
A = 



>)gsx(r-2) 
P 2-r 

MAX 



\ ^max 



- 1 



for r / 2 
for T = 2 



ln C''MAx/ %In) 

Here, Esn it the supernova explosion energy, rj is the frac- 
tion of this energy converted into CRs which are released 
in the interstellar medium, and Emin and Em ax represent 
the extension in energy of the CR spectrum. 

To describe the diffusion of CRs along field lines we 
adopt a diffusion coefficient which is a power-law in energy: 

D *W = B * (lolev) 3 < 8 > 

where D\\ and s are considered here as free parameters. 

Before proceeding in computing the spatial distribution 
of CRs around a SNR we notice that, for z -C Lu, Eq.[2]does 
not provide a good description of the field line wandering. 
The reason is that in this regime the lateral displacement 
of a field line after a step z along Bo is of the order of 
« (SB/Bq)z = bz, since b = SB/Bq represents the angle 
betw een the unpertu rbed (Bo) and total (Bo + SB) magnetic 
field (Isichcnko 1991). An accurate and quantitative analysis 
of the behavior of a magnetic flux tu be in this regim e goes 
beyond the scope of this paper (see llsichenkol Il99ll . for a 
more detailed discussion) . However, in order to describe this 
regime in a qualitative way, for z <C L\\ we substitute Eq. [2] 
with: 



P(R\z) 



[(bz) 2 -(R- Ro) 
%(bz) 2 



(9) 



where $[s] is the Heaviside function, equal to 1 for s > 
and for s < 0. Eq.[9]roughly mimics the behavior of a flux 
tube characterized by a opening angle b. In the intermediate 
region z ~ Lu we use an interpolating function to bridge the 
behaviors described by Eqns. [2] and [9] 
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Figure 1. Cosmic ray over-density around a typical super- 
nova remnant (see text for details) for a particle energy of 
E = 1 TeV at a time t = 10 kyr after the explosion. The left 
panel refers to an isotropic diffusion coefficient of cosmic rays 
equal to D = 5 X 10 26 (£/10 GeV) ' 5 cm 2 /s, while the right 
panel refers to an anisotropic diffusion scenario with Dn = 
10 28 (E/10 GeV) - 5 cm 2 /s, D m = 1 pc, and b 2 = {6B/B ) 2 = 0.2. 
The black cross marks the a position at which the CR over-density 
is equal in the two panels. 



3 RESULTS 

In this Section we compute the spatial distribution of CRs 
expected in the vicinity of a SNR at a given time after the 
explosion. We consider a typical supernova, characterized 
by the following fiducial values of parameters: an explosion 
energy of Esn = 10 erg, a mass of the ejecta equal to 
M e j = 1.4 Mq, and a density of the circumstellar medium 
no = lcm -3 . We further assume that a fraction rj = 0.1 
of the supernova explosion energy is converted into CRs, 
which are injected in the interstellar medium with a power 
law differential energy spectrum dN/dE oc E~ a which ex- 
tends from -Emin = 1 GeV to -Emax = 5 PeV (approximately 
the position of the knee in the CR spectrum). It is known 
from CR data that a shou l d be in the range ~ 2.1 — 2.4 
l|Castellina fc Donatol [20121 : Istrong et ahl 120071 '). We adopt 
here a — 2.2 as a representative value. As described in 
Sec. [5] CRs are gradually released from the SNR during the 
Sedov phase that goes from t w 280 yr to t w 3.6 x 10 4 yr 
ijCioffi et al.lll988T ). For the parallel diffusion coefficient of 
CRs (Eq. HI we assume D\\ = 10 28 cm 2 /s and s = 0.5. 

As a first step, we compare in Fig.[T]the results that are 
obtain ed if an isotropic diffusion coefficient is assume d (as, 
e.g., in lAharonian fc Atovanlll996l : iGabici et alj|2009h . with 
the ones obtained for the anisotropic diffusion model that 
we consider in this paper. In both panels of Fig.[TJ the SNR 
is located at the centre of the field and the color code refers 
to the excess of CRs with respect t o the average density 
of CR s in the Galaxy, which is (e.g. IParticle Data Group I 
l2008h : 



1(E) 



GeV" 



(10) 



Over-densities are plotted for a particle energy of 1 TeV and 
for a time t — 10 kyr after the supernova explosion. Here the 
diffusion coefficient of the magnetic field lines is set equal to 



lpc, with b 1 = (SB/ Bo) 2 = 0.2 (different values of 
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Figure 2. Cosmic ray over— density with respect to the galactic background around a supernova remnant (located in the centre of the 
field). The particle energy is E = 30GeV (upper row), E = 3TeV (middle row) and E = 300 TeV (lower row) and the age of the 
supernova remnant is t = 6kyr (left column), f = 19kyr (middle column) and t = 60kyr (right column), respectively. The diffusion 
coefficients are as in Fig.[T] 



D m will be explored in the following). This corresponds to a 
parallel coherence length of the perturbation of L« — 20 pc. 

The spatial distribution of CRs is strikingly different in 
the two scenarios: spherically symmetric in the left panel, 
and strongly elongated in the direction of the magnetic 
field flux tube in the right panel. A filamentary diffusion 
of CRs was also found in the numerical simulations by 



iGiacinti et ail (|2012T ). The same parameters have been used 
to compute the over-densities in the two scenarios, with 
the exception of the CR diffusion coefficient, which in the 
left panel has been assumed to be isotropic and equal to 
D — D (E/W GeV) 5 cm 2 /s with D = 5 x 10 26 cm 2 /s. The 
choice of two significantly different values for D and D\\, 
with D Du has been made in order to obtain the same 



6 L. Nava and S. Gabici 
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Figure 3. Cosmic ray over— density for particles of 3 TeV around a supernova remnant of age 19 kyr. The assumed values of the 
parameters are as in Fig. [2] except for the diffusion coefficient of the magnetic field lines D m which is 0.5, 1, and 2 pc (left to right 
panel). 



level of CR over-density in the vicinity of the SNR. As an 
example, the black cross in Fig. \T\ shows a position, located 
30 pc away from the centre of the explosion, where the CR 
over-density is identical in the two panels. To get comparable 
values for the CR over-density, a much smaller (isotropic) 
diffusion coefficient D is needed in order to compensate for 
the larger solid angle over which CRs can propagate. As al- 
ready stressed in the introduction, this fact must be taken 
into account when interpreting the gamma-ray observations 
of molecular clouds illuminated by CRs escaping from SNRs. 
This will be discussed in Sec. 0] when the model developed 
here will be applied to fit the gamma-ray observations of 
the SNR W28. 

In Fig.[2]we show the CR over-density around the SNR 
for different values of the particle energy and of the time af- 
ter the supernova explosion. The upper, middle, and lower 
panels refer to a time of 6, 19, and 60 kyr after the explosion, 
respectively. Plots on the first, second, and third column re- 
fer to particle energies of 30 GeV, 3 TeV, and 300 TeV, re- 
spectively. The escape of CRs is described by Eq. [5] which 
states that higher energy CRs are released first, and lower 
energy CRs escape at later times. This is the reason why 
there is no CR excess in the top-left panel of Fig. [2] for 
the choice of parameters made here, for a SNR age of 6 kyr 
particles with an energy of 30 GeV are still confined within 
the SNR shock. As the age of the SNR increases, CRs dif- 
fuse further away along the flux tube and fill a broader and 
broader region. As a consequence of that, the CR over- 
density decreases accordingly. It is evident from these maps 
that a molecular cloud located in the vicinity of the SNR 
can be illuminated by the escaping CRs and become a bright 
gamma-ray source only if it is located within the flux tube. 
A nearby cloud which is not magnetically connected with the 
SNR will not be illuminated by CRs, despite its proximity 
to the SNR. 

All the plots in Fig. [2] refer to a region of size ~ 200 pc 
around the SNR. As said in Sec. [2] this roughly represents 



the e xpected length of a magnetic flux tube in the Galaxy 
(e.g. iPtuskin et all 12008). For distances larger than a few 
hundred parsecs from the SNR, the flux tube loses its iden- 
tity and it is disrupted due to the exponential divergence of 
field lines. Thus, the results presented in this paper are ac- 
curate and reliable for distances up to a few hundred parsecs 
and less. 

Since we assumed here that CRs are strongly magne- 
tized, i.e. their remain attached to field lines, their transport 
across the mean magnetic field is solely governed by the dif- 
fusion of field lines, which is described by the diffusion coef- 
ficient D, n . This quantity determines how quickly the field 
lines diverge as a function of the displacement z along the 
mean field. This is demonstrated in Fig. |3j where the CR 
over-density for particles of energy 3 TeV is shown. The age 
of the SNR is 19 kyr. Three different values for D m are con- 
sidered: 0.5, 1, and 2pc for the left, middle, and right panel, 
respectively. Unfortunately, our knowledge of the properties 
of the interstellar magnetic field is not good enough to allow 
a determination or an estimate of this parameter. For a SNR 
located in a diffuse interstellar medium (i.e. with no massive 
molecular clouds) the morphology of the resulting gamma- 
ray emission due to CR proton-proton interaction with the 
ambient gas would closely follow the spatial distribution of 
CRs. Thus, observing the diffuse gamma-ray emission gen- 
erated by runaway CRs around SNRs might serve as a tool 
to explore the structure of the interstellar magnetic field. 
The detection of such diffuse emission is within the capabili- 
ties of future gam ma-ray instruments such as the Cherenk ov 
Telescope Array l|Acero et al.ll2012l ; ICasanova et alj|2010h . 

Finally, we show in Fig. [4] the spectra of escaping CRs 
at different distances from the SNR and at different times 
after the supernova explosion. Each panel refers to a differ- 
ent epoch: 6, 19, and 60 kyr for the top, middle, and bottom 
panel, respectively. Solid curves show the spectra at three 
different positions on the z-axis: 40 pc, 100 pc, and 200 pc. 
For each of these positions we also show the spectra at differ- 
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ent distances from the z-axis: 25 pc (dotted lines) and 50 pc 
(dashed lines). 

At high energies, in almost all cases the energy spectra 
are power laws with slope w a + s/2, where a is the slope 
of the injection spectrum of runaway CRs and s is the slope 
of the energy dependent diffusion coefficient (Eq. [HJ . Such a 
behavior can be inferred from Eq. [3] If one moves to larger 
values of z, the on-axis (i.e. R = 0) high energy spectrum 
preserves the same slope, but its normalization decreases as 
~ 1/z. This is due to the fact that the transverse section 
of the magnetic flux tube is increasing proportionally to z, 
while the CR intensity along a field line is independent of 
z for z 2 <C 4D||t (see Eqns. [5] and [3}. A feature common 
to all the spectra plotted in Fig. [4] is the presence of a low 
energy cutoff. The cutoff is due to the fact that at a given 
time, only particles of sufficiently large energy had enough 
time to propagate over the distance z. The cutoff is moving 
towards lower energies if a longer time is considered, because 
particles with lower energies have then the time to reach a 
given position z along the axis. Finally, some curves (for 
example the ones with z = 40 pc and R — 25 pc) exhibit a 
quite sharp high energy cutoff. This cutoff is due to the fact 
that particles of different energy are injected within different 
transverse sections of the flux tube. Higher energy particles 
are released earlier from the SNR, when the shock radius 
is smaller, lower energy ones are injected later, when the 
shock radius is larger. While diffusing along the field lines, 
CRs are displaced in the transverse direction due to field 
line wandering. Higher energy particles, which have been 
injected in a smaller region around z = 0, need on average 
a larger transverse displacement in order to reach a given 
distance R from the 2-axis. Thus, for small enough z, the 
opening of the magnetic flux tube might not be enough to 
allow high energy particles to reach R, and this explains the 
presence of the cutoff. 

In the next Section we apply the model developed above 
to a specific object, i.e. the SNR W28 and the molecular 
clouds located in its proximity. Such clouds have been de- 
tected in gamma rays and this has been interpreted by many 
authors as the result of their being illuminated by CRs that 
escaped the SNR. We will demonstrate that a good agree- 
ment can be reached between the predictions of the model 
and observations and we will discuss the impact of this on 
the attempts to derive the particle diffusion coefficient close 
to SNRs by means of gamma-ray observations. 



4 APPLICATION TO THE SUPERNOVA 
REMNANT W28 

W28 is an old SNR in its radiative phase of evolution, lo- 
cated in a region rich of dense molecular gas with average 
density > 5 cm" 3 . At an estimated distance of ~ 2kpc the 
SNR shock radius is ~ 12 pc and its velocity ~ 80km/s (e.g. 



Rho & Borkowski 


|2002| 


Cioffi et al.l ( 1988 


) and 



pernova ejecta is ~ 1.4 Mq, it is possible to infer the super- 
nova explosion energy (Esn ~ 0.4 x 10 erg), initial velocity 
(~ 5500 km/s), and age (t age ~ 4.4 x 10 4 yr). 

Gamma ray emission has been detected from the re - 
gion surrounding W28 both at TeV ( Aharonian et al 120081 ') 
and GeV energies (|Abdo et al.ll2010l ; iGiuliani et alj|20ld) . 



Time = 6x10 3 yr - 




Figure 4. Spectra of runaway cosmic rays at different positions 
and times after the explosion. The age of the supernova remnant is 
t = 6kyr (upper panel), t = 19kyr (middle panel) and t = 60kyr 
(lower panel). Solid lines refer to spectra along the z-axis, oriented 
as the mean magnetic field, at three different positions (z = 40 pc, 
z = 100 pc and z = 200 pc). For each distance 2, three different 
values of R (the perpendicular distance from the z-axis) are also 
considered: Ft = pc, R = 25 pc, and R = 50 pc. The black lines 
show the CR background. 



by HESS, FERMI, and AGILE, respectively. The TeV emis- 
sion correlates quite well with the position of three mas- 
sive molecular clouds, one of which is interacting with the 
north-eastern part of the shell (and corresponds to the TeV 
source HESS J1801-233), and the other two being located to 
the south of the SNR (TeV sources HESS J1800-240 A and 
B). The masses of these clouds can be estimated from CO 
measurements and result in « 5, 6, and 4 x 1O 4 M0, respec- 
tively, and their projected distances from the centre of the 
SNR are w 12, 20, and 20 pc, respectively (|Aharonian et all 
2008). The GeV emission roughly mimics the TeV one, ex- 
cept for the fact that no significant emission is detected at 
the position of HESS J1800-240 A. 

The gamma-ray emission from the clouds in the W28 
region has been interpreted by many authors as the re- 
sult of the interaction of CRs that escaped W28 with the 
dense gas in the cl oud ([Fuiita et all 12009]; ILT fc Chenll2010l ; 
iGabici et all 120101 : lOhira et al.l l201ll ; |Yan et all |2012T ).~A11 
these approaches started from the assumption of isotropic 
diffusion of CRs, and a general consensus was found on 
the fact that, in order to fit observations, the diffusion co- 
efficient had to be suppressed by a factor of w 10... 100 
with respect to the average value in the Galaxy, which is 
D ga i « Dp ( E/lOGeV) 6 with Dp sa 10 28 ...10 29 cm 2 /s and 
s w 0.3... 0.7 |Castellina fc Donatcll20l3 ). 

In this section we take a different approach and we ap- 
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Figure 5. Gamma-ray emission from the three molecular clouds 
surrounding the supernova remnant W28. Fermi data are shown 
as open (blue) circles. Filled (black) circles refer to HESS data. 
The dashed lines show the contribution to the gamma— ray emis- 
sion from the cosmic ray galactic background, the long-dashed 
lines show the contribution from cosmic rays that escaped from 
W28 and the solid (red) line is the total. 

ply the model developed in Sec. [2] to estimate the spectrum 
of CRs and derive the 7-ray emission expected from the 
clouds in the W28 region. This approach is radically differ- 
ent from the ones mentioned above because it relies on the 
more physical assumption that the diffusion of CRs is not 
isotropic, but proceeds mainly along the magnetic field lines. 
The results of our modeling are shown in Fig. [5] where 



the gamma-ray data from the three molecular clouds are 
plotted as blue open symbols (data from FERMI) and black 
filled dots (data from HESS). The emission from the sources 
HESS J1801-233, and HESS J1800-240 A and B is plotted in 
the top, middle, and bottom panel, respectively. The black 
dashed lines represent the contribution to the gamma-ray 
emission from the proton-proton interactions of the CRs in 
the galactic background with the inter-cloud gas. The blue 
long-dashed lines represents the contribution to the emission 
from the runaway CRs that escaped from W28. The solid 
red line is the total emi ssion. The gamma-r ay fluxes have 
been computed following iKamae et all ( 2006) with an addi- 
tional multiplicative factor 1.5 to take into account elements 
heavier tha n hydrogen in both cosmic rays and ambient gas 
l|Morilll997l ). 

A good agreement with observations is obtained is a 
parallel diffusion coefficient — 10 28 cm 2 /s with s = 0.5 is 
adopted, together with a diffusion coefficient for field lines 
D m — lpc with b 2 = (SB/ Bo) 2 = 0.2. Moreover, we as- 
sumed that ~ 20% of the total explosion energy has been 
converted into CRs with a spectrum proportional to E~ 2 ' 2 
and extending from 1 GeV to 5 PeV. In order to be illumi- 
nated by the escaping CRs, the three molecular clouds have 
to be located in the proximity of the axis of the magnetic 
flux tube (i.e. the direction of the local mean field). The 
spectra reported in the figure refers to the positions z = 10, 
165, and 35 pc and R = 6.5, 0, and 14 pc (top to bottom 
panel, respectively). 

It has to be noted that, due to the number of parameters 
involved in the model, other sets of parameter values might 
be found that provide an equally satisfactory fit to data. 
This is not surprising, given that several previous modelings 
of this source provided an equally good fit to data by using a 
radically different picture (i.e. isotropic diffusion of CRs) for 
the transport of particles. Moreover, while a quite small nor- 
malization of the (isotropic) diffusion coefficient, roughly of 
the order of D « 5 x 10 26 cm 2 /s had to be adopted in order 
to fit data satisfactorily, in the anisotropic case we obtain 
a good agreement with data for a significantly larger value 
of the (parallel) diffusion coefficient of ~ 10 28 cm 2 /s. 
It might be noticed that this number is close to the stan- 
dard values inferred for the diffusion of CRs in the Galaxy. 
Thus, any attempt to constrain the CR diffusion coefficient 
from the observations of gamma-rays from the vicinity of 
SNRs needs to take into account that an intrinsic uncer- 
tainty exists, and it is related to the unknown nature of the 
CR transport in the interstellar medium, and in particular 
to the unknown relative relevance of the transport parallel 
and perpendicular to the magnetic field lines. 



5 CONCLUSIONS 

The details of the transport of CRs in the Galaxy are still 
little understood. Studies of the composition of CRs provide 
us with an estimate of the average confinement time of CRs 
within the Galaxy, which can be translated in to a spatially 
averaged diffusion coefficient f or CRs (e.g. IStrong et alj 
l2007l ; ICastellina fc Donatcll2012l ). Whether the CR diffusion 
coefficient has large spatial variations or it is rather uniform 
throughout the Galaxy is not known, thought a suppression 
of diffusion close to CR sources might be expected due to 
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CR s treaming instability (|Ptuskin et al . 2008; M alkov et al.l 
2012). To this purpose, the detection of gamma-ray emis- 
sion from the vicinity of CR accelerators might be used to 
constrain the CR diffusion coefficient , and thus assess the 
importance of such sup pression (e.g. lAharonian fe Atovanl 
1 19961 : iGabici et aD[20 09). This is because CRs escaping the 
accelerators would produce gamma rays via proton-proton 
interactions with the ambient medium. Both the morphol- 
ogy of the resulting emission and its spectrum would de- 
pend on the functional form (i.e. energy dependence, level 
of anisotropy) of the diffusion coefficient. 

An object that has been extensively investigated in this 
context is the SNR W28. Three massive molecular clouds, 
with total mass in the w 10 5 M Q range, are located in the 
vicinity of the SNR shell and emit gamma rays. This has 
been interpreted as the result of the illumination of the 
clouds by the CRs that escaped the SNR. Several mod- 
els have been proposed to fit these observations, and all of 
them are based on the ass umption that the diffusion of CRs 
proce eds isotropically (e.g. lGiuliani et aTll2010l ; lGabici et al.l 
l20ld . and see Sec. [4] for a complete list of references) . There 
is a general consensus on the fact that the (isotropic) diffu- 
sion coefficient has to be suppressed by a factor of m 10... 100 
with respect to the average Galactic one in order to ex- 
plain the observations. This implies coefficients in the range 
D w 10 26 ...10 27 cm 2 /s. 

In this paper, the assumption of isotropy of diffusion 
has been relaxed, and a more physically motivated situa- 
tion have been investigated, in which CRs propagate mainly 
along the magnetic field lines. We considered here the limit- 
ing scenario in which the diffusion of CRs across field lines 
is very small and thus can be neglected. In such a situation, 
the transverse displacemen t of CRs is uniquely d ue to the 
wandering of the field lines (jjokipii fe Parkerlll969T ). Spectra 
and morphology of the spatial distribution of CRs around 
SNRs have been computed and described. The main feature 
is the elongated, filamentary distribution of CRs, as opposed 
to the spherical distribution found in the case of isotropic 
diffusion. 

In order to fit the gamma-ray data from the W28 region 
within this scenario, one has to assume that the molecular 
clouds in its vicinity are magnetically connected to the SNR 
through a magnetic field flux tube. If this is the case, an 
accurate fit to data can be obtained. Under this assump- 
tion, particles are bound to the flux tube and thus forced 
to propagate along a specific direction. For plausible values 
of the diffusion coefficient of magnetic field lines, in order 
to obtain the correct CR over-density at the location of the 
molecular clouds a large (parallel) diffusion coefficient of the 
order of Du « 10 28 cm 2 /s has to be adopted. 

The fact that a very good agreement has been found 
with data in the two radically different scenarios character- 
ized by isotropic and anisotropic diffusion tells us that more 
data needs to be collected from more SNRs in order to in- 
fer with reasonable confidence the properties of the diffusion 
of particles escaping their accelerators. The diffuse emission 
that these runaway particles would produce in their interac- 
tion with the ambient gas is, even in the absence of very 
massive clouds, within the capabil i ties of the Cherenko v 
Telescope Array (|Acero et al 1 120121 : ICasanova et~aH l2010). 
These observations will provide us with precious informa- 
tions about the properties of the transport of CRs in the 



Galaxy, but also with a direct evidence for the fact that 
SNRs are indeed the accelerators of galactic CRs. 
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